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Abstract 



Spin foam models, loop quantum gravity and group field theory are discussed as quantum 
gravity candidate theories and usually involve a continuous Lie group. We advocate here 
I to consider quantum gravity inspired models with finite groups, firstly as a test bed for the 

full theory and secondly as a class of new lattice theories possibly featuring an analogue 
diffeomorphism symmetry. 

To make these notes accessible to readers outside the quantum gravity community we 
^ I provide an introduction to some essential concepts in the loop quantum gravity, spin foam 

and group field theory approach and point out the many connections to lattice field theory 
and condensed matter systems. 

^ ! 1 Introduction 

CN . Spin foam models [H O [3] arose as one possible quantization method for gravity. These models 

can be seen in the tradition of lattice approaches to quantum gravity [4j, in which gravity 
is formulated as a statistical system on either a fixed or fluctuating lattice. Lattice methods 
are also ubiquitous in condensed matter and (Yang Mills) gauge theories. Indeed spin foam like 
structures appear in the strong coupling expansion of QCD or in the high temperature expansion 
for Ising gauge [5] systems. 
IT^ • In any lattice approach to gravity one has to discuss the significance of the (choice of) 

<^ I lattice for physical predictions. Although this point arises also for other lattice field theories, it 

carries extreme importance for general relativity. For matter fields on a lattice we can interpret 
the lattice itself as providing the background geometry and therefore background space time. 
In contrast, in general relativity, the geometry, and therefore space time itself is a dynamical 
variable, and has to be the outcome instead of one of the ingredients of the theory. This is one 
aspect of 'background independence' for quantum gravity [HIT]. 

One possibility to alleviate dependence of physical results on the choice of lattice, is to turn 
the lattice itself into a dynamical variable and to sum over (a certain class of) lattices. This 
is one motivation for the dynamical triangulation approach [8], causal dynamical triangulations 
[9] and group field theory [TOllIl]. See also [12] for another possibility to introduce a dynamical 
lattice. The class of triangulations to be summed over can be restricted for instance in order to 
implement causality [H [13] or to symmetry-reduce models [H] 

Another possibility is to ask for models which are per se lattice or discretization indepen- 
dent. Such models might arise as fixed points of a Wilsonian renormalization flow and repre- 
sent so-called perfect discretizations [15 ^ 1161 [T7] . For gravity this concept is intertwined with 
the appearance of discrete representation of diffeomorphism symmetry in the lattice models 



^In some cases the resulting models might then be reformulated as (tiling) models on a fixed lattice [14) . 
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fIU\ \W[ [TOl I17j . Such a discrete notion of diffeomorphism symmetry can arise in the form 
of vertex translations [2T]: Discrete geometries can be represented by a triangulation carrying 
geometric data, for instance in Regge calculus [22] the lengths of the edges in the triangulation. 
Vertex translations are then symmetries that act (locally) on the vertices of the triangulations 
by changing the geometric data of the adjacent building blocks. If the symmetry is fully imple- 
mented into the model, then this action should not change the weights in the partition function. 

Such vertex translation symmetries are realized in 3-D gravity, where general relativity is 
a topological theory [23]. This means that the theory has no local degrees of freedom, only 
a topology dependent finite number of global ones. The first order formalism of 3D gravity 
coincides with a 3D version of BF theory [24j. BF theory is a gauge theory and can be 
formulated with a gauge group given by a Lie group. 3D gravity is obtained by choosing SU(2). 
However, one can also choose finite groups [251 ESI [27|. In this case one obtains models used 
in quantum computing [28] as well as models describing topological phases of condensed matter 
systems, so called string net models [29] . Moreover BF systems can be seen as Yang-Mills 
like systems for a special choice of the coupling parameter (corresponding to zero temperature). 
Hence in 3-D we have Yang Mills like theories with the usual (lattice) gauge symmetry as well 
as topological BF theories with the additional translation symmetries. 

BF theory (in any dimension) is a topological theory because it has a large gauge group of 
symmetries, which reduces the physical degrees of freedom to a finite number. It is characterized 
not only by the usual (lattice) gauge symmetry determined by the gauge group but also by so- 
called translational symmetries, based on the 3-cells (i.e. cubes) of the lattice and parametrized 
by Lie algbebra elements (if we work with Lie groups). For 3-D gravity these symmetries can 
be interpreted as being associated to the vertices of the dual lattice (for instance given by 
a triangulation) and are hence termed vertex translations. However for 4-D the symmetries 
are still associated to the 3-cells of the lattice and hence to the edges of the dual lattice or 
triangulation. As in general there are more edges then vertices in a triangulation, we obtain 
much more symmetries then the vertex translations one might look for in 4-D gravity. To obtain 
vertex translations of the dual lattice in 4-D one would rather need a symmetry based on the 
4-cells of the lattice. 

In 4-D gravity is not a topological theory. There is however a formulation due to Plebanski 
|30j that starts with -B-F-theory and imposes so-called simplicity constraints on the variables. 
These break the symmetries of -B-F-theory down to a subgroup that can be interpreted as 
diffeomorphsim symmetry (in the continuum). This Plebanski formulation is also the one used 
in many spin foam models. Here one of the main problems is to implement the simplicity 
constraints [311 [32l [33l [M] into the BF partition function. The status of the symmetries is quite 
unclear in these models, however there are indications |20j that in the discrete the reduction of 
the BF translation symmetries does not leave a sufficiently large group to be interpretable as 
diffeomorphisms. Nevertheless if these models are related to gravity, then there is a possibility, 
that diffeomorphism symmetry as the fundamental symmetry of general relativity, arises as a 
symmetry for large scales. In the above mentioned method of regaining symmetries by coarse 
graining, one can then expect that this diffeomorphism symmetry arises as vertex translation 
symmetry (on the vertices of the dual lattice). 

We want to emphasize that hence in 4-D spin foam models are candidates for a new class 
of lattice models in addition to Yang-Mills like systems and -BF-theory. This new class would 
be characterized by a symmetry group in between those for Yang-Mills theory (usual lattice 
gauge symmetry) and -B-F-symmetry (usual lattice gauge symmetry and translation symmetry 
associated to 3-cells). Here we would look for a symmetry given by the usual lattice gauge 
symmetry and translation symmetries associated to 4-cells. 

As we will show in this work the 4Z) spin foam models, which as gravity models are based 
on S'0(4),50(3, 1) or SU{2), can be easily generalized to finite groups (or more generally ten- 
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sor categories [25j). Although the immediate interpretation as gravity models is lost one can 
nevertheless ask whether translation symmetries are realized and if not how these could be im- 
plemented. This question is much easier to answer for finite groups than in the full gravity 
case. One can therefore see these models as a test bed for the full theory. This also applies for 
renormalization and coarse graining techniques which need to be developed to access the large 
scale limit of spin foams. With finite group models it might be in particular possible to access 
the many-particle (that is many simplices or building blocks in the triangulation) and small 
spin (corresponding to small geometrical size of the building blocks) regimqj. This is in contrast 
to the few particle and large spin (semi-classical) regime [36] which is accessible so far. 

In the emergent gravity approaches ^37j one attempts to construct models, which do not 
necessarily start from gravitational or even geometrical variables, but nevertheless show features, 
typical of gravity. The models proposed here can also be seen as candidates for an emergent 
gravity scenario. A related proposal are the truncated Regge models in [38], in which the 
continuous length variables of Regge calculus are replaced by values in Zg, with q = 2 for Tsing 
quantum gravity'. 

Apart from providing a wealth of test models for quantum gravity researchers another intention 
of this work is to give an introduction to some of the spin foam concepts and ideas to researchers 
outside quantum gravity. Indeed spin foams arise as graphical tools in high temperature expan- 
sions of lattice theories, or more generally in the construction of dual models. We will review 
these ideas in section [2] where we will discuss Ising like systems and introduce spin nets as a 
graphical tool for the high temperature expansion. Next, in section[3]we will discuss lattice gauge 
theories with (Abelian) finite groups. Here spin foams arise in the high temperature expansion. 
The zero temperature limit of these theories gives topological BF theories which we discuss 
in more detail in section HI In section [5] we detail spin foam models with non-Abelian finite 
groups. Following a strategy from quantum gravity we will also discuss so-called constrained 
models which arise from BF models by implementing simplicity constraints , here in the form of 
edge projectors, into the partition function. This will in general change the topological character 
of the BF models to non-topological ones. 

We will then discuss in section [6] a canonical description for lattice gauge theories and show 
that for BF theories the transfer operators are given by projectors onto the so called physical 
Hilbert space. These are known as stabilizer conditions in quantum computing and in the 
description of string nets. Here we will discuss the possibility to alter the projectors in order to 
break some subset of the translation symmetries of the BF theories. 

Next, in section [7] we give an overview of group field theories for finite groups. Group 
field theories are then quantum field theories on these finite groups and generate spin foams as 
Feynman diagrams. Finally, in section [8] we discuss the possibility of non-local spinfoams. We 
will end with an outlook in section [9l 

2 Ising like models 

In this section we will review the construction of dual models [391 SO] for Ising like models. 
Similar techniques will lead to the spin foam representation for (Ising) gauge theories. The 
main ingredient is to perform a Fourier expansion of the couplings, which can be understood as 
functions of group variables. Here the group in question is Zq, with q = 2 for the proper Ising 
models. 

^The small spin regime arises as the spin is just a label for representation of the group and for finite groups 
there is only a finite number of irreducible unitary representations. 
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Consider any regular or random lattic^ in which the edges are oriented (the latter is not 
necessary if we consider Z2, that is Ising models). On such a lattice one can for instance define 
an Ising model with spins gy associated to the vertices and nearest neighbour interaction by 
defining a partition function 

9v e 

where = ±1 and group product is the normal product (or addition mod 2 if 5 = 0, 1 is 
used) and /3 is a coupling constant, proportional to the inverse temperature. The action S = 
P 9s{e)9t{e) describes the coupling between the spin g(^(g) at the source or starting vertex and 
the spin (7j(e) at the target or end vertex of every edge. 

More generally we can assume that gy takes any values between 0, . . . ,q — 1 (giving Vector 
Potts models). The group product corresponds now to addition modulo q. A further general- 
ization of (j2.ip is to allow for any edge weigth We (even locally varying) as a function of the two 
group elements gs{e)9~[(^e) associated to an edge e. 

^ = En^^(^^^(e)^^4) • (2.2) 

gv e 

Now the functions We on the group can be expanded in a basis given by the irreducible 
representations of the group [H]. For the groups Zq this is just the usual discrete Fourier 
transform 

q-l 

w{g) = ^w{k)xk{9) 

k=0 

9-1 

w{k) = q-^Y.'^ignM (2.3) 
9=0 

where Xk{9) = exp(^/c • g) are the group characters. Characters for Abelian groups are multi- 
plicative i.e. Xk{9i ■ 92) = Xk{9i) ■ Xk{92) and also Xk{9~'') = Xk^{9) = Xk{9)- 
Applying the character expansion to the weights in (j2.2p we obtain 



J2m2'^e{k)xk{g) 

gv e k 

gv ke \ e / e 



Xkei9) 



= EE n*^(^^) nnx'^e(5°(^'^^) (2.4) 

gv ke \ e / V eDv 

where o{v,e) = 1 if u is the source of e and —1 if v is the target of e. Here we just resorted 

the product such that in the last line the factor Y[eDvXke{9v^^'^^) involves all the appearences 
of gy. We can now sum over gy which according to the Fourier inversion theorem gives a delta 
function involving the of the edges adjacent to v: 

^=5»''E(n*^(^^)) n'^^'^(E'^(^'^)^^) (2.5) 

fc, \ e / v eDv 



^The models can be extended to oriented graphs. We will however not be interested in the most general 
situation but assume that the lattice or more generally cell complex is sufficiently nice in order to construct the 
dual lattices or complexes. Later on we will also need to be able to identify higher dimensional cells (2— cells and 
3— cells) in the lattice. 
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where denotes the number of vertices and 6^'^\-) is the g— periodic delta function. 

We have now a representation of the partition function in which the group variables are 
replaced by variables taking values in the representation labels (that is the Pontrjagin dual of 
G). For Zg the Pontrjagin dual and the group itself are isomorphic, as both g,k = 0, . . . ,q — 1. 
We will call a representation of the form (|2.5p a spin foarcQ representation: there the sum over 
group elements is replaced by a sum over representation labels. Representations meeting at 
vertices (for gauge theories it will be edges) have to satisfy a certain condition, namely that the 
oriented sum of the representation labels has to vanish. This means that the tensor product of 
the representations meeting at an edge has to be trivial. For non-Abelian groups this condition 
is generalized to include some additional data: namely a projector from the tensor product of 
all representations meeting at a vertex into the trivial representation. 

Note that the number of representation labels ke which are associated to the edges is in 
general different from the number of group variables attached to the vertices. But in addition 
we have the restriction that the sum of the representation labels meeting at a vertex has to be 
zero (mod q). These are the Gauss constraints, demanding that the analogue of the electric field 
- namely the assignment of representation labels to the edges, is divergence free. 

To finalize the construction of dual models [JO] one solves the GauB constraints and replaces 
the representation labels associated to the edges of the lattice by variables defined on the dual 
lattice. In two dimensions a divergence free vector field can be constructed from a scalar field 
{vi = d24>, V2 = —di(j)). On the lattice a similar construction leads to a dual model with variables 
associated to the vertices of the dual latticeH Indeed for Z2 one finds the Ising model again, 
just that the high (low) temperature regime of the dual model is mapped to the low (high) 
temperature regime of the original model. 

In three dimensions a divergence free vector field can be constructed from another vector 
field by taking its curl: Vi = tijkdjWj. However, adding a gradient dicf) of a scalar field to 
the vector field Wi will not change the values of the divergence free vector field fj. This will 
be translated into a local gauge symmetry in the dual lattice. For the dual model one hence 
obtains as variables labels on the edges of the dual lattice and the weights are functions of the 
holonomies around the plaquettes (more generally faces, that is the two-dimensional cells) of 
the dual lattice, that is the couplings are between the variables associated to the edges around 
each face. 

This results in a dual model which is a gauge model with variables associated to the edges 
of the dual lattice and with local gauge symmetries defined on the vertices of the dual lattice 
(as in the continuum this gauge symmetry is parametrized by a scalar field). We will discuss 
such theories in more detail in section [3l 

In four dimension, with a similar argument as in three dimensions, the dual model is a 'higher 
gauge theory' ~ the variables are now associated to the plaquettes or faces of the dual lattice 
and the weights define couplings between the plaquettes bounding the three-dimensional cells. 

Finally in one dimension every vertex is adjacent to two edges, hence the Gauss constraints 
in ()2.5p force the representation labels to be equal (we assume periodic boundary conditions). 

^Note that in this particular example we are abusing language twice as we do have neither 'spin' nor a 'foam'. 
The term spin arises in the full models from using the representation labels (spin quantum numbers) of SU (2). The 
proper 'foams' arise for lattice gauge theories as a higher dimensional generalization of the spin nets introduced 
in this section. 

^Here we assume trivial cohomology of the lattice, otherwise one may find so-called topological modes, see 

m 
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We can easily find 

k 

where in the last line we assumed that the couplings are homogenous on the lattice, that is do 
not depend on the position or orientation of the edge. 

The spin foam representation is adapted to the high temperature expansion 03]: Indeed if 
we consider the Ising model with standard couplings as in (12. ip we obtain for the coefficients 
■wqjWi in the character expansion 

wo = cosh(3 , zDi = cosh /3 tanh /3 . (2.7) 

One can now expand the partition function (j2.5|) in powers of wi/wq, that is in the number of 
times the representation label kg = 1 appears. 

Because of the GauB constraints in (12. Sp the number of 'excited' edges, that is carrying 
representation label fcg = 1, entering each vertex has to be even, so for a cubical lattice the 
lowest order contribution is described by one plaquette ~ (wi/wq)'^. This has to be weighted by 
a combinatorical factor, namely the number of ways one can embed the square into the lattice. 

In general the terms in the perturbation series can be represented by closed polymers or 
graphs embedded into the lattice. The expansion in terms of these graphs can then be organized 
in different ways [43]. Consider one particular contribution to the expansion, i.e. a configuration 
where a given number of edges carries non-trivial representation labels. In general such a 
contribution decomposes into connected components of edges with non-trivial representation 
labels, where we define two set of edges with ke ^ as disconnected if these sets do not share a 
vertex. 

We will call such a connected component a (restricted) spin net, where restricted refers to 
the condition, that all edges have to carry a non-trivial representation. Because of the delta- 
functions in ()2.5p such a spin net does not have open ends, i.e. all vertices have either valence 
two or higher. (For the Ising model only even valency is allowed.) Furthermore edges meeting 
at a two-valent vertex have to carry the same representation label. Note that this also holds for 
two-valent vertices, at which the edges 'change direction'. (Here we assume that the couplings 
Ce or We for the edges do not depend on the position or on the direction of the edges. Otherwise 
one has to equip the spin nets with more embedding data.) The representation labels can only 
change at so-called branching points where three or more edges meet. This allows to define 
so-called geometric (restricted) spin netsU which are just specified by the connectivity of the 
spin net vertices (i.e. vertices of valence higher than two with repect to the spin net) and the 
length of the edges of the spin net, that is the number of lattice edges that make up a spin net 
edge. The contribution of such a spin net to the partition function (or the free energy which 
would involve a sum over connected components) can be split into two parts, one is specific to 
the spin net itself, and from (|2.4p given by 

n Mke) , (2.8) 

the other contributions are combinatorical factors determined by the number of embeddings of 
the spin net into the given lattice. These factors could be summarized as measure (or entropic) 

®This definition is taken over from [ij, in which a similar concept was defined for spin foams, that is branched 
surfaces, which we wiU discuss in the next section. The following discussion is motivated by similar considerations 
for 'spin foams' appearing in lattice gauge theories [451 [T] [21 [M] . 
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factors, which carry the information of a given lattice, including its dimension. The spin net 
specific contribution factorizes over the edges of the spin net, for an edge with representation k 
it is given by w{ky where / is the length of the spin net edge (i.e. the number of lattice edges 
that make up this spin net edge). 

The combinatorical factors are determined by the 'background' lattice and its structure. For 
quantum gravity models one would rather be interested in 'background independent' models, 
where such factors should not play a role. Here one can define an alternative partition sum over 
geometric spin nets, whose edges carry representation labels as well as another label, specifying 
the lengths of the edges. Again, in gravity the geometry of space time is a dynamical variable, 
i.e. should rather be encoded in the labels. This motivates to look for models, for which the 
spin net weights are independent of the lengths of the spin net edges so that the weights depend 
only on the labels. These can be termed abstract spin nets [Hj. (Note that with the definition 
we have taken here, the abstract spin nets will in general still remember some features of the 
lattice, for instance the valency of the spin net cannot be higher than the valency of the lattice.) 
For the Ising model such abstract spin nets arise in the limit tanh /9 = 1 for zero temperature. 
Indeed the high temperature expansion is an expansion in the sum of the length of the excited 
edges, whereas the abstract spin networks arise in a regime where this length does not play a 
role. 

A similar structure arises in string net models [29] , which were designed to describe condensed 
phases of (scale free) branching strings. Indeed these string net models are closely related to (the 
canonical formulation of) topological field theories, such as iJF-theories, which we are going to 
describe in section [6l In this canonical formulation spin networks will appear as one choice of 
basis for the Hilbert space of i?F-theories. String net models assign amplitudes (corresponding 
to the so-called physical wave functions in BF-theoiy) to spin networks that satisfy certain 
(gauge) symmetry properties. These imply that these spin networks can be freely deformed on 
the lattice without changing their amplitude. 

In the next section we will consider gauge systems, for which we will obtain a spin foam 
picture. Spin foams are similar to spin nets, just that the role of the edges is taking over by 
plaquettes or faces and the role of the vertices by edges or branchings. Also there one can define 
geometric and abstract spin foams 01] , which are branched surfaces [451 E 12] whose elementary 
surfaces, or spin foam faces, carry a non-trivial representation label. (For non-Abelian groups 
the branchings or spin foam edges carry intertwiners between the representations associated to 
the surfaces meeting at the edges.) The amplitudes for the geometric spin foams will in general 
depend on the area of its elementary surfaces, i.e. the number of plaquettes making up the 
elementary surface in question, and the lengths of the branchings. This is for instance the case 
for Yang-Mills like theories, where geometric spin foams arise in the strong coupling or high 
temperature expansion. 

The conditions for having abstract spin foams, that is amplitudes independent of the area 
of the surfaces, can be understood as requiring independence of the spin foam amplitude under 
trivial (face and edge) subdivision. This can be used to specify measure factors for the spin 
foam models [H 146] . Abstract spin foams will be also generated by group field theories 
discussed in section [71 however the difference is that for the spin foams generated in this way 
trivial representation labels are still allowed. (Also there is not a restriction that these spin 
foams are embeddable into a given lattice.) See also [37] for one possible connection between 
partition functions with restricted (i.e. trivial representations are not allowed) and unrestricted 
spin foams. 

The discussion on geometric and abstract spin nets and spin foams assumes that the couplings 
do not depend on the position or direction of the edges or plaquettes in the lattice, otherwise one 
would need to introduce more decorations for the spin nets and foams. However one can show 
[48j that for instance the Z2 Ising gauge model in AD with varying couplings is universal, i.e. 
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can encode all other %p lattice models. Here it would be interesting to develop the equivalent 
spin foam picture and to see how a spin foam with additional decorations could encode other 
spin foam models. 

3 Gauge theories 

The spin foam representation originated as a representation of partition functions for gauge 
theories [U [H [6] . Such gauge theories can also be formulated with discrete groups [5l HHl [26] 
and for Z2 give the Ising gauge models [5]. A spin foam representation for Yang-Mills theories 
with continuous Lie groups has been presented in |5lj, see also [52], and the results can be easily 
adapted to discrete groups. For abelian (discrete and continuous) groups these are again closely 
related to the construction of dual models [lO] and the high temperature or strong coupling 
expansion |45j . 

Given a lattice (more precisely an orientable two-complex, in which we can uniquely specify 
the 2~cells, called faces or plaquettes) we associate group variables to the (oriented) edges of 
this lattice. We want to construct models with gauge symmetries at the vertices, where as usual 
the gauge action is given by ge — t- gs{e)9eg^(l^ where are gauge parameters associated to the 
vertices of the lattice (remember that s{e) is the source vertex of e and t{e) the target vertex.). 
Gauge invariant quantities are given by (traces in some matrix representation of the group over) 
closed Wilson loops. Expressing the action as a function of these, the action will be gauge 
invariant, as can be easily checked. The 'smallest' Wilson loops are the loops made up of edges 
around a face hf =: Ylge where here we take the oriented product, i.e. the orientation of the 
edges is the one induced by the orientation of the face (and ge = g~}i where is the edge with 
opposite orientation to e). 

Given a unitary (finite dimensional) matrix representation of a group g 1— t- U{g) a (Wilson-) 
action can be defined as 

S = aY,Re(tr{{U{hf))) . (3.1) 

with a a coupling constant. For the groups Zg the irreducible unitary representations are 
one-dimensional and labelled by k = 0, . . . ,q — 1 so that the representations are given by 
g I—)- exp(^ k ■ g). (Here 5 = 0, ... (7 — 1 and the group product is given by addition.) 
More generally we can assume the partition function to be of the form 

^ = iEn^/(^/) (3-2) 

^ 9e f 

where the weights Wf are class functions, i.e. are invariant under conjugation Wf{ghg^^) = 
Wf{h). 

3.1 Spin foam representation 

The (exponential of the) action, and more generally the weights Wf are class function of the 
plaquette variables hf and can hence be expanded in characters. For the Abelian groups Zq this 
again just amounts to the discrete Fourier transform. We can therefore consider the general 
partition function (where we include a measure normalization factor 1/q for every integration 
over the group Zq). 

^ = EinE*^xfc(M • (3-3) 

9e ^ f k 



8 



The derivation of the spin foam representation proceeds similar to the previous section [2J 



9e ^ kf f 

= iE n*fc/ HE n^^/(5' 

kf \ f /eg, /3e 



= E(n-^.)nn^^^^(E°(/'^)^/) (3-4) 

kf \ f J V eDv /De 

where o(/, e) = 1 if the orientation of the edge coincides with the one induced by the (adjacent) 
face, and —1 otherwise. In the last step we just splitted the delta functions over every edge into 
two delta functions over every vertex, as in spin foams one usually works with vertex amplitudes 
(which here would be = HeDu ^^^\'l2fDe ^)^f)- (These amplitudes and the splitting are 
more complicated for non-Abelian groups, see section [5j) 

3.2 Dual models 

Again the construction of the dual models [40] proceeds by solving for the delta functions 
associated to the edges in (j3.4p . that is by solving for the GauB constraints. These enforce that 
the oriented sum of representation labels associated to the faces around one edge is vanishing. 
For the lowest dimensional case, namely two dimensions we obtain a 'trivial' dual theory. Since 
every edge is bounded by two faces, all the representation labels coincide kf = k and we obtain 
(assuming the do not depend on the position and orientation of the plaquette, Wkf = Wk) 

Z = Y,{w,ff . (3.5) 

k 

For the three dimensional case we have seen that the models without any gauge symmetry 
are dual to a gauge model. Hence the Ising gauge model in 31) is dual to the Ising model in 3D. 
In AD the gauge models map again to gauge models. For nontrivial cohomology of the lattice 
topological modes in the dual model might appear |42] . 



3.3 High temperature expansion and spin foams 

As for the Ising models we can also discuss the high temperature expansion, that is for the 
Ising gauge model an expansion in powers of wi/wq = tanh/3. The discussion is very similar to 
the one for the Ising models, the difference being that the perturbative contributions are now 
represented by closed surfaces and not closed polymers any more. 

The 'ground state' is the configuration kf = for all faces. The configurations with some 
kf ^ have to satisfy the GauB constraints, which demands that the number of excited pla- 
quettes (with kf ^ 0) around every edge has to be even. This means in particular that the 
plaquettes with kf ^ have to make up closed surfaces. For instance for the Ising gauge model 
on a hypercubical lattice one would obtain that the lowest order contribution would be repre- 
sented by a {3D) cube giving an order (wi/wq)^ for its six faces (times a combinatorical factor 
giving the numbers of embeddings of the cube into the lattice). 
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More generally the perturbative contributions are now described by closed and branched 
surfaces. (For the free energy one has only to consider connected components t^-) Here a 
proper branching would be a sequence of edges where at least three excited plaquettes meet. 
Excited plaquettes make up the elementary surfaces or spin foam faces. That is at the inner 
edges of these spin foam faces there are always only two excited plaquettes meeting. Again 
because of the GauB constraints, the representation labels of all the plaquettes in a given spin 
foam face have to agree. For homogeneous and isotropic couplings the spin foam amplitude 
(ignoring the combinatorical embedding factors) depends only on the number of plaquettes a 
inside each spin foam face, the contribution of a spin foam face with label k being ~ {wk/wo)'^, 
but not on how the spin foam is embedded into the lattice. (In general for models with non- 
Abelian groups the amplitudes might also depend on the length of the branching or spin foam 
edges.) The definition of geometric and abstract spin foams [44j parallels the one for spin nets. 
Note that the condition for abstract spin foams invalidates the regime of the high temperature 
expansion, which is an expansion in the number of excited plaquettes of the underlying lattice. 

In section m we will discuss BF-theory, which is a topological model and which satisfies the 
conditions for an abstract spin foam, i.e. the amplitude does not depend on the areas a of the 
spin foam faces. The reason is that Wk/wQ = 1 for the SF-models. For the Ising gauge model 
this is the case at zero temperatur^. 



3.4 Expansion around a topological sector 

If we start from the first line in ()3.3p again and keep the summation over both group and 
representation labels we obtain 

9e ^ f k 

= ^ 5^ 2^ n ^M^kf ■ hfige)) 

^ ge kf f ^ 

= :J^Y.Y.]l'''^'P^—''f-^f^9e) + lnvjk^) . (3.6) 

^ 9e kf f ^ 

This corresponds to a first order representation of Yang-Mills theory where ge are the connection 
variables and kf represent the dual (electric field or flux) variables. The model where all = 1 
is a topological BF-model, whose partition function can be solved exactly. Hence Yang-Mills 
theory can be understood as a deformed BF-theoiy [541 [52] - the deformation appears here as 
the ln(wkj:) term (in the continuum it is S A -kB). 

The expansion of models with propagating degrees of freedom, for instance 4Z) gravity and 
Yang-Mills, around topological theories has been discussed for instance in \55\ |3^J 

For the Ising case we can obtain an expansion around the SF-partition function by intro- 

^Indeed in the limit of zero temperature the partition function for the Ising gauge models enforces all plaquette 
holonomies to be trivial. This coincides with the partition function of BF theories discussed in section [H 
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ducing the expansion parameter a = tanh /3 — 1 which is smah for fow temperature: 

z = (cosh/3)tt^^n'^^'^(E°(/'^)^/) (nw^/'0)+'^(^/'i)(i+")) 

kf e /3e \ / 

= (cosh/3)tt/^n'^^'^(E'^(/'^)^/) 

kf e /3e 

^ h /K/2 

... + Jf J2 S{kf„l)---6{kf^^.,l)] . (3.7) 
/i<-</«/ ^ 

That is for the coefficient of we have to count for each configuration {kf}f the number of 
excited faces. For the coefficient of we sum over ordered pairs of faces and in this way count 
for each configuration the number of ordered pairs of excited faces. The next is a sum over 
ordered triples of faces and so on. The last coefficient for a^-^ is one for the configuration kf = 1 
and zero for all other configurations. The terms in this expansion can be seen as expectation 
values of observables in a i?F-model, where the observable is the number of ordered n-tuples 
of excited plaquettes. 

4 Topological models 

In the last section we have seen that the Yang-Mills type theories can be understood as a 
deformation of BF-theories, which we will here discuss in more detail. F stands for the curvature 
of a connection, usually taking values in the Lie group, and i? is a Lie algebra valued {D — 2)- 
form. Hence in constructing these models for discrete groups the question arises with what to 
replace these variables B. We have seen that the field B corresponds to the representation labels. 
Examples for topological models with discrete groups are discussed for instance in [25l [56l [27] . 
See also j57j for the relation between Chern-Simons like theories and Si^-like theories for Zq 
in 3D. In this section we will restrict to Abelian groups, in particular Zg. The models for 
non-Abelian groups will be presented in section [5l 

i?F-theory is a topological field theory in which the equations of motion demand that the 
(local) curvature vanishes. We will start from the partition function encoding this requirement. 
In addition to the usual gauge theory symmetry the model enjoys an additional gauge symmetry 
involving the representation labels, which we will discuss later on. To define the model assume 
that we have given a (oriented) two-complex (or for simplicity a lattice) with oriented edges e 
which bound faces (or plaquettes) /. Associate to the edges group elements ge- (With = g~}i, 
where is the edge with reversed orientation to e.) Then we start from the following definition 
for a partition function Z (see for instance [21 [58] ) 

^ = E^ell^oihf) (4.1) 

9e ^ f 

where again hf = ITec/^'e ™6ans the oriented product of edges around a face. Here we normal- 
ized the integration measure over the group to one |l = 1 where q is the cardinality of the 

group. 6g is the delta function with respect to this measure, i.e. Ylg \^G{9)f{g) = /(id). Note 
that this is just a special case of the gauge theory partition functions (j3.2p . by setting Wf = 6g 
there. 
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The partition function can be easily evaluated to 



Z 



configurations {gele) such that T T = id V / 

J-J-eC/ 



J/ 



where the second factor arises due to the normalization factors for the group integrations over 
Q'e and for the delta functions for each face. 

For the (Abelian) cyclic groups Zq we obtain by using the expansion be = X^fc Xfc 

where o(/, e) = 1 if the orientation of e agrees with the one induced from / and —1 otherwise. 
The term in the exponential can be interpreted as a i?F-action, with the i3~field given by the 
assignment / i— )■ A;^ of representation labels to faces and the curvature F by the holonomies 
/ij = riec/^e- (In the exponential we used the representation of as additive group.) 

The spin foam representation can be derived in the same way as for the gauge theories in 
section I3.lt 

= EH n^^'^(E'^(/'«)M ■ (4-4) 

This gives Z as 

Z = (J [configurations {/c J } J, such that o(/, e)A;j = mod g V e] (4.5) 

that is as the number of assignements {kf}f satisfying the GauB constraint for every edge. In 
the next section we will discuss the local symmetries of the BF partition function. We will find 
that on the {kf}f there is a translation symmetry acting which is based on the 3-cells of the 
lattice. Hence Z refiects the orbit volume of this translational gauge symrnetry. This leads to 
divergencies for the BF partition functions for (continuous) Lie groups [59] 



4.1 Symmetries of the partition function 

Here we will discuss the gauge symmetries of the partition function ()4.ip . As it is of the form of 
a gauge theory (13. 2p the partition function is invariant under the usual gauge transformations: 
9e — ^ gs(e)9eg^(ly where G G are gauge parameters associated to the vertices. This is obvious 
in the representation (|4.ip for the partition function. 

^However not all divergencies are due to this gauge symmetry [59] . 
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There is a further symmetry, the so-called translation symmetry, which is easiest to see in 
the spin foam representation. This symmetry is based on the 3— cells c of the lattice, i.e. the 
cubes for a hypercubical lattice^ Consider a field c i— )■ /cc of gauge parameters associated to the 
cubes and define the gauge transformations 

k'f = kf + ^o{c,f)kc modq (4.6) 

where o(c, /) = 1 if the orientation of / agrees with the one induced by c and —1 otherwise. Then 
if {kf}f is a configuration which satisfies all the GauB constraints, so is {^/}/ and vice versa. 
The reason is that if an edge e is in the boundary of a 3— cell c then there are exactly two faces 
/, /' which are both in the boundary of c and adjacent to e. The orientation factors are such, 
that the contributions of k^ from the two faces /, /' are of opposite sign in the GauB constraint 
for e. Hence its value does not change, and therefore the contribution of the two configurations 
{kf}f,{k'j}f to the partition function is equal. 

This is the principle, that 'the boundary of a boundary is zero' dod = 0, which underlines the 
Bianchi identity for the curvature. With the help of this Bianchi identity one can also explain 
the translation symmetry, see \23\ [18] . For 3D gravity this translation symmetry corresponds to 
the diff'eomorphism symmetry of the theory. The gauge field kc corresponds in the dual lattice 
to a field associated to the dual vertices and one can interpret the gauge transformation as a 
translatioiJ^ of these dual vertices (which in the gravity models are the vertices in a triangulation 
of space time). 

In AD the 3— cells are dual to edges in the dual lattice and the translation symmetries rather 
translate these instead of the dual vertices. However for gravity-like models one would need to 
break down the translation symmetries from the BF-case, where it is based on the dual edges, 
to symmetries based on the dual vertices. (Correspondingly in 4D diffeomorphism symmetry 
of the action leads to Noether charges encoding the contracted Bianchi identities and not the 
Bianchi identities itself.) In the case of gravity one strategy is to impose so-called simplicity 
constraints [30l [311 [32l [Mj into the partition function. For the models with finite groups the 
question arises, if there is a class of AD gauge models with (translation like) symmetries based 
on the dual vertices. Such models would be non-topological and feature propagating degrees of 
freedom, as the symmetry group based on dual vertices is much smaller than for BF, where it 
is based on dual edges. See also the discussion in sections [5] and [6l 



5 Gauge theories for non— Abelian groups 

In the following, we will consider how to generalize the construction performed so far for finite, 
but non- Abelian groups G. Details about representations can be found in pT|. 

Again, we denote by p the irreducible, unitary representations of G on C" ~ Vp where 
n = dim/?. Note that for non- Abelian groups n can be larger than 1. Every function / : G — )• C 
can be decomposed into matrix elements of representations, i.e. 

fia) = ^V^i^fpabP{9)ab (5.1) 
p 

where the sum ranges over all (equivalence classes) of irreducible representations of G. The fp 

2D this symmetry degenerates into a global symmetry: the Gaufi constraints force all kf to be equal, 
kf = k, however the contribution of every representation label k to the partition function is the same. 

^''the gauge parameters are then lie algebra valued fields and for the lie algebra sii(2) indeed correspond to 
translations 
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are given by 



fp,a,b = j^^J2V(^^^P fia) p{9)ab (5.2) 

where |G| denotes the number of elements in the group G. Introducing an inner product between 
functions /i , /2 : G — )■ C by 

(/l|/2> = TLYjf^(9)f2(.9) = J2^^^^P'<'^f'{f2)p,a,b (5.3) 

' ' gee p,a,b 
then the matrix elements of p are orthogonal, i.e. 

{p.,\Pm) = . (5.4) 

Furthermore, we denote by 

Xp{g) = tT{p{g)) (5.5) 
the character of p, then the J-function on the group G is given by 

^cig) = ^dim pxp{g) (5.6) 
p 

which satisfies 

i^^^SGig)fig) = /(I) (5.7) 
5.1 G— gauge theory on a two-complex 

Consider an oriented two-comple K. Denote the set of edges E and the set of faces F, then 

by a connection we mean an assignment — )• G of group elements ge to edges e £ E. For every 
face f £ F, we denote the curvature of this connection by 

hf ■■= II 9°^^'^^ (5.8) 

eedf 

the ordered product of group elements of edges in the boundary of /, where we take o(/, e) = ±1 
depending on the relative orientation of e and /. In the non-Abelian case the ordering of the 
group elements around a face/plaquette actually matters, and we choose here the convention as 
depicted in figure [H 

As before we define a gauge invariant partition function by choosing a collection of weight- 
functions Wf : G ^ C invariant under conjugation. These encode the action (and path integral 
measure) of the system. The partition function is defined as 

' ' 9e feF 



^^More specifically, we consider only complexes which are such that the glueing maps used to successively build 
up K axe injective. This in particular means that the boundary of each face contains each edge at most once. 
With certain choices for amplitudes, this condition can be relaxed slightly, see e.g. ^S] 
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ei 



Figure 1: A face /, bordered by edges ei, . . . , 65. The curvature is given by /i/ = ye4S'e3 ye2S'ei 
(note the ordering). 



where |G| is the number of elements in G. Since wj can be expanded into characters via (|5.ip . 
the path integral can be written as 

' ' 9e Pf f 

In the sum there is for each edge e a representation Pf{ge) appearing for every face / adjacent 
to e, / D e. The corresponding sum results in (assuming the orientations of all / agree with 
that of e for the moment, in order not to overburden the notation) 

{Pe)hi2-i„; jlj2-jn ■= T7r^^Pfli9)hjlPf2i9)i2j2 Pf„i9)inj„ (^-H) 

' ' geG 

It is not hard to see that the operator P, called the Haar-intertwiner, given by 

{Pe'^)i\i2---in '■— {-Pe)iii2---in; jlj2---jn'^jlj2---jn • (^-12) 

which maps the edge-Hilbert space 

Tie ■■= Vp,C^Vp^(g) ...C^Vp„ (5.13) 

to itself, is actually an orthogonal projector on the gauge-invariant subspace of He- El Note 
that, while in the Abelian case one has to sum over all representations over faces such that at 
all edges the (oriented) sum adds up to zero (as in (|3.4|) ). in the non- Abelian generalization 
one has to sum over all representations such that at each edge e the representation of the faces 
/ D e meeting at e have to contain, in their tensor product, the trivial representation. Also, one 
obtains a nontrivial tensor Pg for each edge, which in the case of Abelian gauge groups is just 
the C-number (5±A:i±fc2± - ±fc„,o- Since the representations for non-Abelian groups can be more 
than 1-dimensional, in general the tensor Pg has indices which are contracted at each vertex v 
in the two-complex k. 

Choosing an orthonormal base Li^\ A; = 1, . . . , m for the invariant subspace of Tie, we get 



k=l 



In case the orientations of e and / do not agree for some e, then is essentially replaced by 
in ()5.1ip . which leads to the appearance of the dual^ representation in the tensor product 



^^This can be easily shown by proving Pe = , resulting from the unitarity of all representations, and (Pe)^ 
Pe, which uses the translation-invariance and normalization of the Haar measure on G. 
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Defined by /9*(/i)a6 = ^)ba- 
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(j5.13p of the He- In this case be labels a basis of intertwiner maps between the tensor product 
of all representations associated to faces / with o{f,e) = — 1 and the tensor product of all 
representations associated to the faces with o{f,e) = 1. 

In (j5.14p . one regards |z.e) and (te| as attached to, respectively, the endpoint and the beginning 
of e. For each vertex v in k, the associated |te) and (te| can be contracted in a canonical way: 
For every face / which touches v, there are exactly tho edges ei, 62 in the boundary of / that 
meet at v. The definitions above are exactly such that, if one chooses a basis in each Vp and the 
dual basis in each V*, in and the two indices associated to / are in opposite position, so 
the can be contracted. See figure [2] for an example. 




Figure 2: The relative orientations of ei and /, and 62 and / agree, so in both the edge Hilbert 
spaces Tiei and 7ie2 the factor Vp^. Therefore |iei) and (te2l have opposite indices belonging to 
Vp,, which therefore can be contracted. 



Therefore, contracting all the appropriate Le at one vertex leaves one with the vertex- 
amplitude Av{pf, Le), which depends on the representations pf and intertwiners Le associated to 
the faces and edges that meet at v (and the orientations of these) . The vertex amplitude can be 
computed by evaluating the so-called neighbouring spin network function, living on graph which 
results from a dimensional reduction of the neighbourhood of v: Construct a spin network, where 
there is a vertex for each edge e touching v, and a line between any two vertices for each face / 
between two edges. Assigning the Pf to the lines of the spin network, and the intertwiners |te) 
or {Le\, to the vertices, depending on whether the corresponding edge e is incoming or outgoing 
of V, results in a spin network, whose evaluation gives Av 

With this, the state sum can, using ()5.10p and ()5.1ip . be written in terms of vertex amplitudes 

via 

^ = E n(*/)p ll^v{Pf,Le) (5.15) 

Pf l-e f V 

5.2 Examples 

The first example we consider is G-BF theory, which corresponds to an integral over only 
flat connections, i.e. the choice Wf{h) = Scih). Since the 5-function can be decomposed into 
characters with {wf)p = dim pf, we get 

Z = J2 Il^^^^Pf ll-^v . (5.16) 

pf Le f V 

If K is the two-complex dual to a triangulation of a 3-dimensional manifold, the vertex amplitude 
is essentially the analogue of the 6j-symbol for G0. 

^*Not that there is some ambiguity in the Uterature about what is actually called the 6j-symbol, which is a 
question of the correct normalization. We mean the normalized 6j-symbol here, since we chose the intertwiners 
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The second example one usually considers is Yang-Mills theory. The Wilson action can be 
specified, similar to the Abelian case, by choosing a unitary representation /), so that 

Svuih) = I {Xp{h) + Xp{h-^)) = «Re(xp(M) (5-17) 
where a is the coupling constant. The weights are then given by Wf{h) = exp(— S'yM (/l))• 
5.3 Constrained models 

There is a generalization of the state-sum models (|5.15|) . coming from the desire to obtain 
non-top ological generalizations of BF-theoiy. It amounts to choosing the intertwiners Lg not to 
span all of the invariant subspace, but only a proper subspace Vg C Inv(^e)- This originates 
in the attempt to define a state-sum model for general relativity, which can be written as a 
constrained BF-theoiy. The subspace Ve is viewed as the space of intertwiners satisfying the 
so-called simplicity constraints, see e.g. [301 [601 [M] . Examples for such models are the Barrett- 
Crane model [31] and the EPRL and FK models [32]. These models can also be written in the 
form of a path integral over connections [6H [STl [62] , but we shall not concern ourselves with 
this here. 

In the following, we will introduce a class of models which can be seen as the generalization 
of the Barrett-Crane models to finite groups. Given a finite group H, the model is a state-sum 
model for the group G = H x H. Irreducible representations of G are then pairs of irreducible 
representations {p^ ,p^) of H. In the edge Hilbert-space He there is a specific element lbc of 
the space of gauge-invariant elements, called Barrett- Crane-intertwiner. It is only nonzero iff 
the representations and are dual to each other. In that case, for an edge e with attached 
faces /fc, consider the Hilbert space 

Ue = (5.18) 

Again, we have assumed that the orientations of all faces fk and e agree, in order not to 
overburden the notation, otherwise replace by its dual, or equivalently, exchange p^ . If we 

consider the projector Pg : %(. — )• %(, onto the subspace of elements invariant under the action of 
H, then can be seen as an element of Tie ®'H*ei which is naturally isomorphic to Tie-, because 
Pj are dual to each other. The corresponding element lbc in T~^e is invariant under H x H, as 
one can readily see, and it is our distinguished element, onto which Pg is projecting. 

Since the projection space for each edge is (at most) 1-dimensional, the vertex amplitude 
for the BC-model does not depend on intertwiners, but only on the representations p^ = (p/)* 
associated to the faces. In fact, since the representations are unitary, every representation is 
dual to itself, so we effectively have only one representation pj attached to each face. Using 
diagrammatical calculus [63j, it is not hard to show that the vertex amplitude for the BC model 
is given by a sum over squares of the BF-theoiy model on the same vertex, i.e. 

= Y.\-^rHpf,^e)\' (5.19) 

where the sum ranges over an orthonormal basis tg of S'3-intertwiners in Tie for every edge e at 
the vertex v. 

Let us shortly discuss the case oi H = Ss, the group of permutations in three elements. For 
53, which is generated by (12), the permutation of the first two elements, and (123), which is 

to be normalized in the first place. The normalization, which usually results in nontrivial edge-amplitudes, can 
be absorbed into the vertex amplitude. Also there might, according to convention, some sign factors be assigned 
to the edges e, which may not be absorbed into the vertex amplitudes Av 
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the cyclic permutation, subject to the relation (12)^ = (123)^ = 1, the representations theory is 
well-known ^41j: There are three irreducible representations, called [1], [—1] and [2]. The first 
is the trivial representation, the second one maps a permutation a to its sign (—1)°". The third 
one is two-dimensional, and 

"(t^^)) - ( ; -1 ) ■ ""^^^o - ( 2is "ci ) f^-^") 

The nontrivial tensor products of the representations decompose as 

[-1]0[-1] = [1], [-1]»[2] = [2] (5.21) 
[2]0[2] = [l]e[-l]©[2] 

Therefore the intertwiner space in Tie is e.g. three-dimensional, when there are four faces at- 
tached to e carrying the representation [2]. Hence, for H = S3, the BC-model is an example for 
a constrained version oi a H x iif-state-sum model, as defined in the last section, because Pg 
projects onto a proper subspace of the invariant subspace of Tie- 

This class of models is an example for an abstract spin foam [44], since its amplitudes only 
depend on combinatorial data, and the topology of the two-complex. In particular, it leads to a 
model which is invariant under refinement of the two-complex k, by trivial subdivisions of edges 
or faces. Hence these models provide potential examples for models which are background- 
independent, but not topological (as is the case in the full theory). 

Note that there is a wealth of different models, which come from different choices of nontrivial 
subspaces of He, onto which Pe in projecting. In particular, one could consider EPRL-like models 
where G = S4 x S4 or G = A4 x A4 and consider the subspace of intertwiners for p'^ = pj, 
which are in the image of a boosting map b : V^^ — )• Vp+ V~. , given by the fusion coefficients 
(see e.g. [32j for H = SU{2)). Here ^4 (^4) is the group of (even) permutation of four elements. 
Another, possibly more geometric way would be to consider embeddings ^4 — ^5 and thus 
construct proper subspaces of Ss-intertwinersEl Such models can be considered as truncations 
of the EPRL models as (^4) is the (chiral) symmetry group of the tetrahedron and a subgroup 
of the rotation group. Here it will be interesting to investigate the relation to the full models in 
more detail. 

The models can serve to test many proposals for the full theory, for instance how the different 
choices of edge projectors Pe determine the physical degrees of freedom or particle content of 
the given theory. This is related to the problem of implementing the simplicity constraints into 
spin foam models. We are planning to investigate the features of these models, in particular its 
symmetries |64] and behaviour under coarse graining, in future work. 

6 Hilbert space, transfer operators and constraints 

In this section we are going to derive the transfer operators [651 [50] for the Yang-Mills like gauge 
theories with partition function ()3.2|5.9p with a finite group G with cardinality \G\. Here we 
will also allow for non-Abelian groups. The first part of the discussion will be similar to [50] for 
Yang-Mills theory, however we will also include SF-theory as a special case. From the transfer 
operators one can obtain in a limiting procedure the Hamiltonian operators. For iJi^'-theory 
this limiting procedure will not be necessary, as we will see, there the transfer operators can be 
understood as projectors onto the space of so-called physical states, which can be characterized 

^^We are thankful to Frank Hellmann for this insight. 
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as being in the kernel of (Hamiltonian) constraints. This illustrates the general principle that a 
path integral with local symmetries acts a projector on a constraint subspace |66j. On the other 
hand, one can construct a projector onto the constraint subspace as a path integral. See [67] . 
in which this is performed for 3D gravity, which can be formulated as an SU (2) — BF-theoiy. 

The transfer operator T is defined on a Hilbert space T-L, so that the partition function can 
be written as 

Z = tr^f^ . (6.1) 

Here we assume for simplicitj0 that we have a hyper-cubical lattice, with one lattice direction 
designated as time direction, in which we have periodic boundary conditions. The trace tr-^ is 
the trace in a Hilbert space which we will describe in the following. The Hilbert space is then 
associated to the spatial lattice and can be built as a direct product of Hilbert spaces associated 
to the spatial edges Cg, so that Ti = T^e^- As the (configuration) variables associated to the 
edges are group elements, a Hilbert space associated to an edge will be Tie = €[(7], i.e. the space 
of functions on the group, equipped with an inner product (all functions have finite norm as we 
consider finite groups), which we will give below. 

We will work in the coordinate (that is holonomy or connection) representation, and choose 
as basis states eigenstates \g). For any unitary matrix representation g i— )• Pabid) of the group 
G these are common eigenstates of the so called edge holonomy operators p{ge)ab 

P{9e)ah\9) = P{9e)ab\g) ■ (6.2) 

These basis states are orthonormal 

(5'l<7)=n^«(4'9eJ (6.3) 

where 5g is the delta function on the group such that ^cls'i 5') ~ 1- have also the 

completeness relation 



IdH = |^El5)(5l . (6.4) 

We use this resolution of identity times in (|6.ip to rewrite the partition function as 

1 ^ 

Z = tr^T^ = — ^ J]n^5n+i|T|5n) (6.5) 

ges,nn=0 

where we introduced n = 0, . . . , to label the 'constant time hyper surf aces' in the lattice. On 
the other hand we will assume that our partition function is of the form 

1 1^1 

^ = ^En^/(M = n ^/.(^/^ n^/*(^/*) u^u^fs) .(6.6) 

' ' 9- / ' ' n=0 (/.)„+! (/«)„ (/.)„ 

Here we introduced weight functions Wf of the holonomies around plaquettes hf, which in the 
form on the right hand side can also chosen to differ for spatial plaquettes fg and timelike 
plaquettes ft- (This is necessary if one wants to obtain the Hamiltonian in the limit of small 
lattice constant in timelike direction.) Since we are dealing with a gauge theory the weights Wf 
are class functions, i.e. invariant under conjugation, furthermore we will assume that Wf{h) = 



^^See [68lll9j for a possibility to introduce a slicing of triangulations into hyp ersurf aces, and a transfer operator 
based on so-called tent moves. 
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Wf{h^^), i.e. the weights are independent of the orientation of the face. A weight function 
can therefore be expanded into characters, which are hnear combinations of matrix elements 
of representations (in fact characters are just the trace). Hence the weight functions wih be 
quantized as edge holonomy operators. 

On the right hand side of (j6.6p we have just spht the product into factors labehed by the 
time parameter n. Here we assume that the plaquettes {ft)n are the ones between time n and 
n + 1 . Comparing the two forms of the partition functions ()6.5p and (j6.6p we can conclude that 



{9n+i\T\gn) 



fs 



' ' 9et (ft) 



n 

fs 



(6.7) 



The two outer factors can be easily quantized as multiplication operators, so that we can 
write 



WKW 



(6.8) 
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Figure 3: A time like plaquette in a cubical lattice. The group elements at the time like edges 
can be understood to act as gauge transformations on the vertices either at time steps n or 
n + 1. Integration over the group elements assoicated to the time like edges enforces therefore 
(via group averaging) a projector onto the gauge invariant Hilbert space. 



To tackle the operator K, note that the group elements associated to the timelike edges 
appearing in the plaquette weights 

^/t (5{e.)„5(e0„5(-e;)„^,5("4)„ ) (6.10) 

can be understood as acting as gauge transformations associated to the vertices of the spatial 
lattice, see figure [31 

That is we define operators r(7), 7 G G^""^ by 

f(7)l5) = \l'^>9) (6.11) 

where (7 > g)e = ^s{e)9e%(l^ and s(e),t(e) are the source and target vertex of the edge e, 
respectively. The operators T generate gauge transformations as 

{glthm = {7>9m = H7>9) ■ (6.12) 
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Now we can see the plaquette weight in (j6.10p as either a wave function at time n or a wave 
function at time n + 1 on which the group elements associated to the timelike edges act as a 
gauge transformation. The sum in (j6.9p over the group elements g^t then induces an averaging 
over all gauge transformations, i.e. a projection onto the space of gauge invariant states. Hence 
we can write 

k = PgKo = KoPg (6.13) 

where 

Pg = T^^^il) > ((7n+l|i^0|5n> = Yl^fMe^Ua^,])^,^^) ■ (6-14) 

' ' 7 ft 

The operator Pg is a projector, that is Pq = Pg- One can easily show that T{'y)PG = pGi'i'y) = 
Pg for any 7 E G'^'"^ , hence it projects onto the space of gauge invariant functions. As the 
operator W is made up from gauge invariant plaquette couplings it commutes with Pg, so that 
we can write 

f = PgWKoWPg . (6.15) 

Here we see an example for a general mechanism, namely that the transfer operator for 
a partition function with local gauge symmetries acts as a projector onto the space of gauge 
invariant states. We can characterize such gauge invariant states as being annihilated by con- 
straint operators. In the case of the usual lattice gauge symmetries these constraints are the 
GauB constraints, which we will discuss later-on. 

To discuss the action of Kq we introduce first the spin network basis, in which Kq will be 
diagonal. 



6.1 Spin network basis 

So far we encountered the holonomy operators, i.e. the multplication operators p{ge)ab in the 
configuration representation. The conjugated operators are the electric fields or fluxes, which 
act as (matrix) multiplication operators in the spin network basis. This is a convenient basis to 
discuss the remaining part Kq of the transfer operator. 

To obtain the spin network basis [UlT] we just need to use that every function on the group can 
be expanded as a sum over the matrix elements p{-)ab of all irreducible unitary representations 
p. For Abelian groups all irreducible representations are one-dimensional, hence a,b = and 
we obtain the discrete Fourier transform (j2.3p . In the general case of non- Abelian groups we 
define spin net work states \p,a,b) by 

Jl Vdimpe Pe{9e)a,b, = {g\p,a,b) . (6.16) 

e 

which are orthonormal, i.e. {p, a, b\p' , a', b') = SppiSaa'^bV 

In the spin network basis we can easily express the left Le{h) and right translation operators 
Re{h). These are the conjugated operators to the holonomy operators 16. 2[ To simplify notation 
we will just consider the Hilbert space and states associated to one edge and omit the edge 
subindex. We define 

L{h)\g) = \h-^g) , R(h)\g) = \gh) . (6.17) 
In the spin network basis 

{9\Hh)\p,a,b) = {hg\p,a,b) = ^/ dim p p{hg)ab (6.18) 
= y^dim p p{h)acp{g)cb = p{h)ac{g\p,c,b) . 
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where we sum over repeated indices. That is 

L{h)\p,a,b) = p{h)ac\p,c,b) , R\p,a,b) = \p,a,c) p{h~^)cb ■ (6.19) 

The remaining factor Kq of the transfer operator factorizes over the edges and it is straight- 
forward to check that its acts on one edge as 

(Here one has to use that Wf^ is a class function and invariant under inversion of the argument.) 
On a spin network state we obtain 



Ko\p,a,b) = ■^^'^Wf^p{h)ac\p,c,b) 
1 

]G\ 



1^1 . 

h p' 

where in the second hne we used that a class function can be expanded into characters xp ^-^id 
in the third the orthogonality relation 



1^ J^Xp'(Mp(^)ac = -^^^P,P'^<^^ ■ (6-22) 

h 

The expansion coefficients Cp are given by 

" = |G| 



^p = ^\Y^ ^ft ih)Xp{h) (6.23) 



h 



That is Kq acts diagonal in the spin network basis and as the eigenvalues only depend on 
the representation labels p it commutes with left and right translations. For Yang-Mills theory 
(with Lie groups) in the limit of continuous time one obtains for Kq the Laplacian on the group 
[651 [50]. 



6.2 Gauge invariant spin nets 

Given a spin network function ip, which can be regarded as function ip{he) of holonomies associ- 
ated to edges, where h G G^'^ is a connection, and a gauge transformation 7 G G"^ an assignment 
of group elements to vertices of the network, then the gauge transformed spin network function 
is given by 

f(7)V'(/le) := V'(7s(e)/ie7t7e)) (6-24) 

where s(e) and t{e) are the source and the target vertices of the edge e. A gauge transformation 
can therefore be expressed as a linear operator in terms of L and R, as shown in the last section. 
Decomposing the spin network function into the basis \pe,ae,be), i.e. 

V'(^e) = ^ '>l'p,,a,,bJPe,ae,be) (6.25) 

Pe,a,e,be 
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one can readily see that the condition for a function ■0 to be invariant under all gauge transfor- 
mations ag translates into a condition for the coefficients '4'pe,ae,be: namely 

Vipe,a... = nil...,.e„ (6-26) 

V 

where the product ranges over vertices v, and each tensor which has the indices ae for 
each edge e outgoing and 6(e) for each edge e incoming v, is an invariant element of the tensor 
representation space 

i^"^ e Inv I (g) Vp^ (g) j (6.27) 

\i)=s(e) v=t{e) j 

i.e. an intertwiner between the tensor product of incoming and outgoing representations for each 
vertex. An orthonormal basis on the space of gauge-invariant spin network functions therefore 
corresponds to a choice of orthonormal intertwiner in each tensor product representation space 
for each vertex, where the inner product is the tensor product of the inner products on each V^, 

p 

Note that neither the holonomies p{he)ab nor left- and right translations are gauge-invariant, 
therefore they map gauge-invariant functions to non-gauge-invariant ones. However, they can 
be combined into gauge-invariant combinations, such as the holonomy of a Wilson loop within 
a net. 



6.3 Local constraint operators and physical inner product 

We have seen that for a general gauge partition function of the form (16. 6p the transfer operator 
(I6J5]1 

f = PgWKqWPg (6.28) 

incorporates a projection Pq onto the space of gauge invariant states. This is a general feature 
of partition functions with gauge symmetries [66]. Indeed in quantum gravity one attempts 
to construct a projector onto gauge invariant states by constructing an appropriate partition 
function. 

As has been discussed in section 14.11 the BF~modeh enjoy a further gauge symmetry, the 
translation symmetries (14. 6p . (A generalization of these symmetries holds also for non-Abelian 
groups.) Hence we can expect that another projector will appear in the transfer operator. Indeed 
it will turn out that the transfer operator for the BF-models is just a combination of projectors. 

This is straightforward to see as for the BF-models the plaquette weights are given by 
Wf{h) = dcih), so that 

W = UiScihjjf^ . (6.29) 

Here one might be worried by taking the square roots of the delta functions, however we are on 
a finite group. The operator is diagonal in the basis l^f) and has only two eigenvalues: q"^-^" 
on states where all plaquette holonomies vanish and zero otherwise. The square root of is 
defned by taking the square root of these eigenvalues. 

Hence W is proportional to another projector Pp, which projects onto the states for which 
all the plaquette holonomies are trivial, that is states with zero (local) curvature. 
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The final factor in the transfer operator T is Kq which according to the matrix element of 
in (I6.14P (putting Wf{h) = 5g(^)) is proportional to the identity 



Ko = IGl^"' Idn (6.30) 

where the numerical factor arises due to our normalization convention for the group delta func- 
tions, which amounts to 5G(id) = Hence the transfer operator for i?F-theory is given 
by 

f = PgPf = PfPg (6.31) 

and since for the projectors we have Pq = PcPp = Pp the partition function simplifies to 

Z = tinf^ = IGf^^^'+^^'hinPcpF =■■ \G\^^^^'+^^^^ dimCHphys) ■ (6.32) 

The projection operators in the transfer operator ensure that only the so called physical states 
ip G Tiphys are contributing to the partition function Z. These are states in the image of the 
projectors (we will call the corresponding subspace Tiphys) and can be equivalently described to 
be annihilated by constraint operators, which we will discuss below. 

The objective in canonical approaches to quantum gravity [6l [7] is to characterize the space 
of physical states according to the constraints given by general relativity, which follow from the 
local symmetries of the theory. (The quantization of these constraint in a consistent way is highly 
complicated [69].) Indeed also in general relativity, as well as in any theory which is invariant 
under reparametrizations of the time parameter, one would expect that the transfer operator 
(or the path integral) is given by a projector, so that the property ~ T would hold. 0rhis 
would also imply a certain notion of discretization independence, namely, the independence 
of transition amplitudes on the number of time steps used in the discretization, see also the 
discussion in [17j on the relation between reparametrization symmetry in the time parameter 
and discretization independence. For a lattice based on triangulations one can introduce a local 
notion of evolution, the so called tent moves [68^ I19j. These moves evolve just one vertex and 
the adjacent cells. Local transfer operators can be associated to these tent moves. These would 
evolve the spatial hypersurface not globally, but locally, and would allow to choose some arbitrary 
order of the vertices to evolve. In this way one can generate different slicings of a triangulation 
as well as different triangulations. The condition that the transfer operator associated to a tent 
move is given by a projector would then implement an even stronger notion of triangulation 
independence. 

However reparametrization invariance, which would lead to such transfer operators, is usually 
broken by discretizations already for one-dimensional toy models. For (coarse graining and 
renormalization) methods to regain this symmetry, see \17\ \70\ I16j . In the case of gravity, 
or more generally non~topological field theories, these methods will however lead to non-local 
couplings for the partition functions |70] . for which one would have to modify the transfer 
operator formalism. 

Furthermore one has to construct a physical inner product on the space of the physical sates. 
This process is quite trivial in the case considered here, as we are considering finite groups and 
all the projectors are proper projectors. Hence the physical states are proper states in the (so- 
called) kinematical Hilbert space Ti and the inner product of this space can be taken over for the 
physical states. In contrast, already for BF-iheory with (compact) Lie-groups the translation 
symmetries lead to non-compact orbits. This leads to physical states which are not normalizable 

^^In the limit of taking the discretization in time direction to the continuum we would obtain tlie same projector 
as for finite time discretization. Indeed the projectors aheady encode for finite time steps the (Hamiltonian) 
constraints, which are the generators of time translations (time reparametrization symmetry). 
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with respect to the inner product of the kinematical Hilbert space. For the intricacies of this case 
(with SU{2)) see for instance [6^. A further comphcation for gravity is the very comphcated 
non-commutative structure of the constraints |71| . 

Let us summarize the projectors into P = Pq Pf- Also in the case of generalized projectors 
a physical inner product can be defined 1 72\ [7] between equivalence classes of kinematical states 
labelled by representatives ipi,'4>2 through 

ii^Mphys := . (6.33) 

Kinematical states which project to the same (physical) state define the same equivalence class. 
This leads for instance to an identification of the two states (for the group Z2) in figure HI This 
is analogous to the action of the spatial diffeomorphism group in (3-D and AD) loop quantum 
gravity, see the discussion in [67J and refiects the concept of abstract spin foams, whose amplitude 
does not depend on the embedding into the lattice, in section [2l Indeed any two states which 
can be deformed into each other by applying the so-called stabilizer operators introduced below 
in (|6.34l6.40p are equivalent to each other. The reason is that the physical Hilbert space is the 
common eigenspace to the eigenvalue 1 with respect to all the stabilizer operators. Hence any 
part of a state that undergoes a non-trivial change if a stabilizer is applied, will be projected 
out in the physical inner product. 
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Figure 4: Two states in the spin network basis for Z2 which are equivalent under the physical 
inner product. The thick edges carry non-trivial representations. The right state can be obtained 
from the left state by applying a plaquette stabilizer. 

Another problem in quantum gravity is then to find (quantum) observables [7S] that are 
well defined on the physical Hilbert space. In the case of gravity these Dirac observables are 
very hard to obtain explicitly, even classically there are almost none known . In particular, 
such Dirac observables have to be non-local [Tl]. If one interprets a quantum gravity model 
as a statistical system a related task is to find a well-defined order parameter characterizing 
the phases. Expectation values of gauge variant observables may be projected to zero under the 
physical inner product [75]: As in our case we work with finite dimensional Hilbert spaces and the 
physical Hilbert space is just a subspace of the kinemtaical one, there is a construction principle 
for physical observables. For any operator O on the kinematical Hilbert space one can consider 
POP, which preserves the physical Hilbert space and corresponds to a 'gauge averaging' of O. 
However many observables constructed in this way will turn out to be constants. For i?i^-theory 
observables can be constructed by considering a product of the stabilizer operators discussed 
below along non-contractable loops \28\ 176]. The resulting observables are non-local, and the 
cardinality of an set of independent operators (equivalently the dimension of the physical Hilbert 
space) just depends on the topology of the lattice and not on its size. 

We will now discuss the stabilizer operators which characterize the physical states. In topo- 
logical quantum computing these are known as stabilizer conditions [28]. 

We have considered the operators r(7) in (16. lip that implement a gauge transformation 
according to the assingements (7)1,^ of group elements to the vertices of the (spatial sub-) 
lattice. These can be easily localized to the vertices Vg of the spatial lattice by choosing the 
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gauge parameters 7 (an assignement of one group element to every vertex in the spatial lattice) 
such, that all group elements are trivial except for one vertex Vg- We will call the corresponding 
operators r^,^ (7) where now 7 denotes an element in the gauge group G (and not in the direct 
product Gl*^=). These can be expressed by the left and right translation operators (j6.17p 

f..(7)b) = n n ^e'Mia) ■ (6.34) 



Physical states lip) have to satisfj0 



f..(7)IV'> = IV') (6.35) 



for all vertices Vg and all group elements 7. As the operators T define a representation of the 
group we can restrict to the elements of a generating set of the group. 

This suggest to consider the action of these (star) operators in the spin network basis: 

^Vs{l)\P,a,b) = Yi /'e(7)aeCe n Pe'{l'^)d^,b^, \P,{Ce,ae',ae"),{be,de',be")) (6.36) 
e:s{e)=Vs e':t(e')=Vs 

where on the right hand side edges e are the ones starting at Vs, e' those ending at Vs and e" 
are all the remaining edges in the spatial lattice. From this expression one can conclude that 
the spin network states transform at Vg in a tensor product of representations 

Pv, = (g) Pe €D (g) pl, (6.37) 

e:s{e)=Vs e':t{e')=Vs 

where p* is the dual representation to p. Gauge invariant states can be constructed by considering 
the projections of these representations to the trivial one, or equivalently by contracting the 
matrix indices of the states with the intertwiners between the tensor product of 'incoming' 
representation and the tensor product of 'outcoming' representations at the vertices, see section 



For the Abelian groups Zq the irreducible representations are one-dimenionsional, hence we 
can omit the indices a, b in the spin network basis. The tensor product of representations (j6.37p 
is also one dimensional and equal to the trivial representation, provided 

^ fee = Yl (6-3^) 

e:s{e)=V3 e':t(e')=Vs 

for the representation labels k^, /cg/ of the outgoing and incoming edges at Vg respectively. Hence 
we recover the Gaufi constraints from the spin foam representation (j3.4p . Here these Gaufi con- 
straints appear as based on vertices as opposed to edges as in (13. 4p . But the spatial vertices Vg 
represent the time like edges and the GauB constraints in the canonical formalism are just the 
Gaufi constraints associated to the time like edges in the spin foam representation (|3.4p . 

Let us turn to the other projector Pp. It maps to states for which all the spatial plaquettes 
fg have trivial holonomies. Hence the conditions on physical states can be written as 



^*Here f „^ is a unitary operator, so that the condition on the physical states is in the form of a so-called 
stabilizer. Alternatively the condition can be expressed by self-adjoint constraints Cu, (7) — i(r„^(7) — ri,^(7~^)) 
for group elements 7 with 7 7^ 7~^. Otherwise define (7) = (7) — Idn- Physical states are then annihilated 
by the constraints. 
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where we have introduced the plaquette operators (corresponding to the flatness constraints) 



Here p should be a faithful representation, otherwise one might find, that also states with local 
non-trivial holonomies satisfy the conditions ()6.39p . see for instance the discussion in |82]. On 
the other hand this can be taken as one possible generalization of the model. For non-Abelian 
groups the plaquette operators additionally depend on the choice of a vertex adjacent to the 
face, at which the holonomy around the face starts and ends. However the conditions (j6.39p 
do not depend on this choice of vertex: if a holonomy around a face is trivial for one choice of 
vertex it will be trivial for all other vertices in this face, as these holonomies just differ by a 
conjugation. 

i?F-theory in any dimension is a topological field theory, that is there are only finitely many 
physical degrees of freedom depending on the topology of space. Consequently the number of 
physical states, or equivalently of equivalence classes of kinematical states in the sense ()6.33p is 
finite and does not scale with the lattice volume. 

There are a number of generalizations one could consider. First one can couple matter in the 
form of defects to the models. In 2)D SU{2)-BF theory corresponds to a first order formulation 
of 2>D gravity. Point-like particles can be coupled to the model, see for instance [76] and lead to 
the violation of the fiatness constrainst ()6.40p (coupling to the mass of the particles) and to the 
violation of the Gaufi constraints (I6.35P (coupling to the spin of the particles) at the position 
of the particle. The defects can be understood as changing the topology of the spatial lattice, 
i.e. to change its first fundamental group. To have the same effect in, say four dimensions, one 
needs to couple strings, see [77] . 

In the context of quantum computing [28] one does not necessarily impose the conditions 
(j6.35p and (|6.40p as constraints but as characterizations for the ground states of the system. 
Elementary excitations or quasi particles are states in which either the flatness or the GauB con- 
straints are violated. These excitations appear in pairs (at least for Abelian models [78j) and 
can be created by so-called ribbon operators [281 ESI ES] > which in the context of the proper 
particle models in 2>D gravity are related to gauge invariant (Dirac) observables |80] . For further 
generalizations based on symmetry breaking from G to a subgroup of G, see for instance [78]. 

In AD gravity is not topological, rather there are propagating degrees of freedom. Similar 
to the discussion for the partition functions in 14.11 one can try to construct new models which 
are nearer to gravity by breaking down the symmetries of i?F-theory. In 3D the flatness 
constraints are defined on the plaquettes of the lattice. Constraints act also as generators of 
gauge transformations (indeed the conditions (I6.35|6.40j) impose that physical states have to be 
gauge invariant) and the fiatness constraints can be interpreted as translating the vertices of 
the dual lattice in space time [81 ^ 118 ^ [20]. which can be seen as an action of a diffeomorphism. 
In AD the same interpretation holds only for some exceptionally cases corresponding to special 
lattices that do not lead to space time curvature, see the discussion in ^60j. In general the 
fiatness constraints rather generate translations of the edges of the dual lattice [83] . 

A possible generalization leading to gravity-like models is therefore to replace the fiatness 
conditions (|6.39p based on plaquettes with some contractions of these conditions, such that 
these new conditions are based on the 3-cells, i.e. cubes in a hyper-cubical lattice. This 
would correspond to the contractions of the form FEE (Hamiltonian constraints) and EE 
(diffeomorphism constraints) in the (complex) Ashtekar variables of the curvature E with the 
electric fiux variables E [84[ [6l [7] . The difficulty in constructing such models is consistency of the 
constraint algebra, i.e. to find constraints that form a closed algebra. However to consider such 
models for finite groups should be much simpler than in the case of full gravity. Even if such 




(6.40) 
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consistent constraint algebras cannot be found, physical Hilbert spaces could be constructed, 
for instance with the techniques in \7T\ [85] . 

Furthermore it will be illuminating to derive the transfer operators for the constrained mod- 
els introduced in section 15.31 as in the full theory the relation between the covariant models and 
Hamiltonians in the canonical formalism is still open [86]. To this end a connection representa- 
tion [62] can be employed. 

7 Group field theory (GFT) 

We have seen in earlier sections that the spin foam formalism is the natural description of the 
high temperature expansion for random lattices. This is by no means the only instance that 
random lattices have been used to analyze the properties of dynamical systems. As we shall see, 
group field theories [11] describe weighted sums over n-dimensional random lattices and have 
enough freedom to capture the dynamics of many dynamical systems. They attempt to capture 
the degrees of freedom of local spin foam models stemming from a sum over cellular complexes. 

In the 2-dimensional case, they are known by the more familiar moniker of (random) matrix 
models |87[ [88] . These are well-studied in a variety of fields and can describe a multitude of 
different scenarios including, but not limited to, 2d gravity with cosmological constant [89], 2d 
Ising/Potts models |90j . 2d Yang-Mills theories |91] . certain string theories ^92], the enumeration 
of virtual knots and links [93], the list goes on. Moreover, they have inspired the development 
of a plethora of useful techniques to solve and analyse statistical ensembles of random matrices 
[M] such as the topological expansion [95] [96], the eigenvalue method [97], the double-scaling 
limit, the method of orthogonal polynomials [98] and the character expansion [99], to name just 
a few. The more recent label of GFT is to higher-dimensional random lattice theories as matrix 
model is to 2d random lattice theories and they share the same motivation. The ambition is to 
develop a similar technology for GFTs and the development of an apparatus to count the degree 
of divergence [lOlj and later a topological expansion has had its first success [102] . but tools to 
extensively solve these models are still missing. 

To say it in yet another way, GFTs are a natural tensor generalization of matrix models. 
They generate a weighted sum over n-dimensional cell complexes in the perturbative expansion of 
the free energy. They have garnered increasing attention from the quantum gravity community, 
with the advent of the spin foam formalism. Not only do spin foam diagrams naturally arise as 
the Feynman graphs of GFTs [1] but inherent in all current spin foam models for 4d quantum 
gravity is a truncation of the degrees of freedom and a manifest loss of diffeomorphism symmetry. 
It is with the hope of re-obtaining these lost properties (e.g. perhaps in a the continuum limit 
regime) that one analyses these higher-dimensional theories. Tensor models first appeared in the 
quantum gravity literature with the work of Ambjorn et al. [100] . Later they were generalized to 
genuine multi-variable field theories over compact Lie groups, such as SU(2) (or more generally 
the internal local symmetry group of the n-dimensional gravity theory in question). These are 
known as Boulatov-Ooguri field theories [10[ I103[ 11041 1105] . With this further generalization, 
one does not think of these models as generating discrete cell complexes with unit n-volume 
cells. Rather the field arguments carry the local dynamical geometrical information for the 
cell complex and thus each Feynman graph already contains a restricted sum over geometries, 
at least in theory. The sum over graphs completes the sum over geometries (and topologies). 
Moreover, higher-dimensional finite group models deserve more study since we shall see that 
they are instrumental in describing certain matter degrees of freedom. 

To commence, let us detail some of the general set-up. 
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7.1 General set-up 



While one might hope to proceed at the outset with the most generic definition possible, this is 
probably not the wisest course of action. On the other hand, we shall not be too specific either, 
but try only to put forth just a flavour of the theory at this stage. With this in mind, let us 
define a somewhat, but not totally, general GFT |106] . The fundamental dynamical object is 
the field (p ■ G^" — )• R; a field on n copies of a group G. Schematically, the action takes the 
form: 



G 

A 

n 



Jl dga dcja 
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(7.1) 



v(K.%i}) 



We use a normalized measure on G, while fC and V are the kinetic and vertex operators. As usual, 
they are instrumental in the construction of Feynman amplitudes. One should pay attention 
to the peculiar coupling of the arguments in ()7.ip . especially in the vertex operator. As a field 
theory on G^", the vertex interaction is highly non-local. But therein lies the power of the 
GFT framework. When one calculates the perturbative expansion of the free energy, one can in 
fact view the Feynman graphs as n-dimensional cell complexes F. What is more, for the 0""''^ 
interaction above, this cell-complex is topologically dual to an n-dimensional triangulation A. 

To get an idea of how this might be the case, let us start with a generic Feynman diagram 
occurring in the expansion of the free energy: 



F = \ogZ 



V(j) e 



-5M 



1 



\Aum\ 



Zr where Zr = X*"" J /xr(G) JJ P JJ V (7.2 



eer ver 



and ^r(G) is shorthand for the measure on F, |Ant[r]| is the order of the discrete automorphism 
group of the graph and V = K,~^ is the propagator. To explain, let us make use of Fig. [5l 
The diagram F is simply an n + 1-valent graph with components given on the first line of the 
diagram. But due to the non-local nature of the vertex operator, we can define a related strand 
diagram siV) which captures more accurately the richness of the information encoded in the 
theory. (Strand diagrams are also known as a fat or ribbon graphs.) For the strand diagram, 
we draw the field rather as a collection of n strands. The coupling described in V({(?^ ^ 9b a}) 
drawn as a re-routing of the strands within the vertex interaction. The vertex we chose in (j7.ip 
contains enough information to be interpreted as the dual 2-skeleton to an n-simplex. In this 
dual picture the field represents a (re — l)-simplex. Its re arguments are the (n — 2)-simplices on 
the boundary of the (n — l)-simplex. The vertex re-routing identifies (re — l)-simplices along their 
(n — 2)-simplices so as to form an n-simplex. Thus, the interaction may be equivalently seen as 
generating elementary n-simplices. The kinetic term simply glues these n-simplices together by 
identifying a pair of (re — l)-simplices on the respective boundaries. To sum up, the Feynman 
diagrams with their amplitudes Z-p are rich enough in content to deserve a strand diagram, 
which in turn may be viewed as the dual to an n-dimensional simplicial complex. One then has 
an arena within which to frame many of the interesting questions of lattice theories. 
Obviously, there are a number of ways that one can generalize the formalism above: 

• It is in fact rather rare that one chooses a real map. Often, one takes a complex field 
^ : G^" — )• C , but in principle there are very few restrictions. In fact, it is convenient to 
view the field as carrying a representation of some other group H (finite or compact). We 
shall see that this allows one to incorporate additional degrees of freedom, interpretable 
as matter. 
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Figure 5: For the case n = 3, we have drawn the fundamental components for the Feynman 
diagram F, the strand diagram s(F) and the dual n-simplex A. 



• One could also include more than one field. This is to a certain extent a sub-case of the 
previous point but a specific type known as colouring has attracted a large amount of 
interest |107j . It has been shown that in the n-dimensional model (17. ip , if one has n + 1 
fields : G^" — )• C, {i G {1, . . . ,n + 1}), and a specific type of interaction Lagrangian: 



one has much more control over the topology of the discrete structures generated in the 
perturbative expansion. For n = 2, one is guaranteed to generate triangulations corre- 
sponding to oriented connected Riemann surfaces (although, in this case colouring is not 
required as long as the field is complex). Even with colouring, it is no longer the case in 
higher-dimensional examples that the triangulation corresponds to a manifold. In n = 3 
for example, one gets not just manifolds, but also objects known as pseudo-manifolds. The 
benefit of having a coloured model is that it removes those pseudo-manifolds from the sum 
that have been shown to contribute more heavily than manifolds to the free energy |108j . 

• Lastly, but certainly not leastly, the presence or absence of symmetries plays a prominent 
role in GFTs. Usually, these symmetries come in a rather idiosyncratic fashion. Some 
models maintain a symmetry under the action of the permutation group on the arguments 
of the field. Positing this symmetry serves to increase the number of topologies arising in 
the sum and/or alters the relative weight of existing configurations. 

In matrix models, anti-symmetry or hermiticity (of the matrix) is a desirable property from 
the point of view of solvability via the method of eigenvalues. There are higher-dimensional 
versions, which have not been linked as yet to methods of solution. 

By far the most important symmetry in n-dimensional GFT is the so-called diagonal 
symmetry under the action of the group. It is these symmetry generators that often truly 
parameterize the dynamical degrees of freedom of the model at the Feynman diagram level. 
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7.2 Matrix models 



7.2.1 2d gravity with cosmological constant 

Consider an Hermitian N x N matrix: (pg^g^ sucli tiiat {4>^)gag^ = {(F')gagb — ^gagt- make the 
Unk to group field theories more transparent, we shall define a related field: (f) : Z^- x Zat — )• C 
and a operation f on the space of fields such that (p^ = cj) where 4>^{ga,gb) = 4'{9b,9a)- Then, the 
action is usually taken to be: 

■S[(t>] = <Piaa,gb) 4>{9a,gb) - ^ (p{9a,gb)(l>{gb,gc)(t>{gc,ga) 

9a, 9b 9a,gb,9c 

X] (t>{ga,gb)(t>{gb,9c)<P{9c,ga) (7.4) 

9a,9b,9c 

= trN{<p')-^trN{(ty^). 
V iV 



where tr^ is the usual trace over N x N matrices and we have used hermiticity to simplify the 
form of the action. Let us examine the free energy: 

where: 

IVCP= fUdcPg^g^ n dd\c{cl>g^g,) dJmicPg^g,) Z = \* ^ N^^'^'^ (7.6) 

ga ga,gb-ga<gb 

The perturbative expansion of the free energy generates a sum over 2d orientable connected 
triangulations. We shall say more about the Feynmanology in a moment. We are in a quantum 
mechanical setting so we can define the measure precisely by taking the Lebesgue measure on 
each of the matrix components. Since orientable Riemann surfaces are classified by a single 
parameter, the Euler characteristic x(^) = xO^Tl)! we can consider the expansion as a sum 
over Riemann surfaces TZ and tri-valent cellular decompositions T-ji. Thus, the matrix model 
naturally encodes a sum over topologies. We construct the Feynman amplitudes as follows. 
Due to the matrix structure of the fields, the strand diagrams associated to the propagator and 
vertex operator are, see Fig. [6l 

9b 9a 



9a » gn 



9c 9c 

Figure 6: Diagrammatics for (p^ hermitian matrix model. 

It is easy to see that these fat graphs s{Tn) are topologically dual to 2d triangulations At^. The 
vertex operator generates the triangles and the propagator glues them pairwise along edges. 
From here, one can calculate the amplitude associated to a given triangulation in the sum. 
There is a factor A / ^/N for every vertex v G T-ji and for every face of Ft^ there is a redundant 
sum which contributes a factor N. Thus, the contribution of the Ft^ to the free energy is: 

Zr^ = X*" Ar#/-i#^ = X*" iv#^~#e+#/ = x*f atxW = x*^ N^-^sn (7.7) 
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where v, e and / are the vertices, edges and faces of T-ji, respectively; which satisfy 3#v = 2j^e; 
and finahy x(^) is the Euler character which satisfies x(^) = xi^n) = #v — #e + #f = 2 — 2g'ji, 
where gji is the genus of the Riemann surface. 

This amphtude looks familiar, however, from 2d quantum gravity. In two dimensions, the 
Einstein-Hilbert action with Riemannian metric is integrable: 

Snig] = j^d\^{R[g] - 2A) = -ax{n)+l3An (7.8) 

where g, R[g] are the metric and Ricci scalar; G and A are the gravitational and cosmological 
constants; while a = 1/4G and /3 = A/SvrG. Thus, on a fixed manifold, 2d gravity is trivial. In 
the quantum case, however, it is no longer trivial, since large quantum fluctuations may change 
the genus of the surface and on higher genus surfaces, there are non-trivial topological sectors. 
Thus, one should define the partition function for 2d quantum gravity as: 

Z2d = J2 [ ^fe^'^^lfl, where X = Riem[n]/Diff[n] (7.9) 

with Riem[}Z]/ Diff[}Z] the space of Riemannian metrics modulo diffeomorphisms. If one con- 
siders the matrix model to generate triangulations with every triangle assigned unit area, then 
^A-ji = #/• Thus, equating: 

N = e^, A = e-^ / = E i /| ,!r 11 ' ^^'^O) 

Jx jr^ \Aut[rn\\ 

we see that the free energy of the matrix model gives the partition for 2d gravity: 

-^ = ^2d (7.11) 

Furthermore, one can see immediately from the topological expansion that the 2-sphere domi- 
nates in the large- A'^ limit (also called the planar limit for that reason): 

J- ""^"^ Z^^s^ = y ,, ]^ (7.12) 



More interesting is the so-called double scaling limit. It turns out that for large numbers of 
vertices: 

^Zr^, ^'^i^^^^ (A - A.)^-- where ,>2 (7.13) 

Thus, it diverges for some critical coupling and we can extract a continuum limit by tuning 
A — )• Ac- We refer the reader to |87j for details. 

The presence/absence of critical behaviour is one of the most striking and important features 
of any matrix model. In fact, for the one-hermitian matrix model, we can generalize to other 
potentials (eg. n-valent graphs dual to n-gons). The continuum limit still exists and describes 
pure gravity once again [109J. 

Furthermore, this model is an example of a class of matrix models which have been solved 
exactly, in this case, using the hermiticity of the matrix to recast the problem in terms of its 
eigenvalues. As a final point, which is of particular interest when studying quantum gravity 
regimes, one can examine the Schwinger-Dyson equations for this model, which generate a loop 
algebra, which is in turn a part of the Virasoro algebra [llOj . 
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7.2.2 2d Ising and Potts models 



The Ising model has a natural realization as a two-matrix model, where the two matrices repre- 
sent the lb states of an Ising spin system. Once again, we can think of it as a field (p : Znx'^n ^ ^ 
which also carries the defining Z2 representation. Thus, (j) = (ph^, where /ij G Z2. The dynamics 
are specified by the action: 



1 



(7.14) 



Setting c = e , we find that the free energy takes the form: 



and A 



EE 



sym[TTi]\ 



Zr^ where Zr^ = 



exp 



f/'E" 



(7.15) 



Ijz^r^. We see that the perturbative expansion generates 4-valent Feynman graphs 



which are dual to quadrangulations of Riemann surfaces. The action contains two types of 
vertex interaction, representing the two spin states of the system, while the propagators couple 
the vertices in a nearest neighbour fashion. Thus, we see that the amplitudes are that of an 
Ising model on the lattice T-j^. 

It turns out that the Ising model is much easier to solve when summed over random lattices. 
Moreover, it is solvable even in the presence of a non-zero magnetic field, see [90]. This ease 
stems from the extra symmetry gained by coupling to gravity, since the difficulties inherent in 
dealing with a fixed lattice are integrated out. 

This form of matrix model can be extended to the ^-state Potts model rather easily. Now 
take q N X N Hermitian matrices, or equivalently a field (p : Zn x Zn ^ ^ which also carries 



the defining representation of Zq, so that 
generalization of the Ising form: 



where hi G Z„. The action is a simple 



N 



(7.16) 



and 



/ 1 



-Cg \ 



1 ... 



where 



e2/3 + g _ 2 • 



(7.17) 



\ -Cg -Cg ... 1 / 

While one can still examine such features as critical behaviour, these models have not been 
solved exactly although there are arguments to suggest that this is possible [lllj . 



7.3 More general GFTs 

7.3.1 n-dimensional topological models: Boulatov-Ooguri models 

Currently, one possesses a less well developed arsenal with which to attack higher-dimensional 
models. All the same, to start off, we can check to see what classes of spin foam amplitudes 
can be reproduced and to what dynamical systems they can be associated. For the topological 
models set out in Section [H we can proceed directly to the GFT formalism. We consider a field 
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C The action in this case is: 



Y[ dga dcja 



.a=l 



(t>{9l,---,9n) 



.a=l 



n+1 

H dga,b 

a,b=l:bj^a 



'n+1 

n 'l>(9a. 
.a=l 



a-1, 



<t>{9i,- 

9a,a—n ) 



,9n) 



n+1 



n ^(9a,b9bi] 
a,b=l:b^a 



(7.18) 

where unhke matrix models, it is conventional to use a normalized measure on the group (these 
normahzed integrals are replaced by normalized sums in the finite group case). The key ingre- 
dient is that the fields possess a symmetry alluded to earlier: 



4>{9i ,---,9n)= (t>{9i9, 9n9) for any G G. 



(7.19) 



We enforce this symmetry by explicit use of a projector onto the space of invariant fields, which 
results in the following propagator and vertex operator: 



V 



dg 



G 



'[lSi9a99a^] 



.a=l 



V 



QXn + 1 



n+l 
.a=l 



n+1 



n ^i9a,b 9a^ 9b 9f,,a) 
a,b=l:bj^a 



(7.20) 



By construction, the amplitude for a given Feynman graph: 

Zr,BF = A*" y n HHGf) where = J] 5° 



(7.21) 



eer 



eedf 



Thus, we see that we have generated the amplitude associated to an n-dimensional BF theory 
()4.3p . This class of models is to some extent the starting point for the analysis of higher- 
dimensional GFTs. For G either a finite or Lie group, the Feynman amplitudes have been 
well studied and their properties documented \102\ I119j . Various aspects have been analysed, 
such as the insertion of topological matter defects [112j, fermionic matter coupling |113j . the 
first instantonic solutions have been found and analysed, a dual representation has been 
proposed |116j along with a first class of radiative corrections [115| . Since these theories are 
background independent, it is natural to expect some trace of diffeomorphism symmetry jTTT] 
and the first steps in the implementation of mean field theory techniques have been completed 

mi- 

Most outstandingly, it appears we can still capture, in these tensor models, certain aspects of 
the topological expansion we saw for matrix models. Indeed, at leading order, for the coloured 
models [107j, the n-sphere dominates \102\ 1119] , Although, these are all interesting topics and 
it would be wonderful to dwell on them in more detail, our aim is to give a more basic outline 
of how GFTs can capture interesting dynamical systems. With this in mind, let us pass to the 
higher-dimensional nearest neighbour models. 



7.3.2 n-dimensional Ising/Potts models and more general nearest-neighbour mod- 
els 

There seems to be a rather general algorithm to generate Ising/Potts models on higher-dimensional 
random lattices. This is simply to copy the 2d case. Let us consider q complex fields cph '■ G^" — )• 
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C where /i € Zg and with the action as before: 



W dga dga 



a=l 



QXn(n + l) 



n+1 
a,b=l:b^a 



\_a=l 

n+1 

4'hi (,9a,a—li • • • ) 9a, a— n 



fPhAgi, ■■■,9n) 



.a=l 



n+1 
a,b=l:b^a 



(7.22) 

In this case, one gets a Potts model coupled to the BF theory outlined in the previous section: 



BF 



n E 



exp (^Y.^<^A{e)j^ 



(7.23) 



where the coupling constant occurring in (|7.23p is the A of equation (j7.15p . In fact, we can 
generalize to even more general nearest-neighbour amplitudes just by altering the propagator 
K^i^o , Consider H a finite (or compact) group. Essentially, choosing: 

K'^^''^ = {we{h.^f)y^ (7.24) 



is enough to generate amplitudes of the form: 



Zv = Zy^bf 



HE 



eer 



(7.25) 



Of course, for these nearest-neighbour interactions, the matter degrees of freedom are always 
coupled to some other theory, such as a BF theory in this case. This is exactly analogous to 
the 2d scenario. 



7.3.3 n-dimensional Yang-Mills and more gauge theory models 



A different approach must be taken for the models of Section [3l It is more tricky to generalize 
to amplitudes of the form: 



Zr = / JJ dge JJ ^fiGf) where Gf 
eer /er 



e&df 



(7.26) 



These are very much in the i?F-type (j7.2ip . but they need not be topological. As mentioned 
before, they may be viewed as a modification of the topological amplitude. We have seen that 
lattice Yang-Mills amplitudes are in this form. But the GFT formalism is not at first sight the 
most natural setting for these amplitudes. This stems from the fact that we do not have control 
over face amplitudes at the outset. The faces / G F arise after gluing of the interaction vertices. 
We do have control, however, over the amplitude assigned to part of the face. This part is known 
as a wedge in the spin foam literature |120j and we shall denote it /„. Generically, the best we 
can do is to specify the vertex operator to be: 



Then, the amplitudes turn to be: 



'n+1 



.a=l 



n+1 



n ^hi9a,b9a^ 9i,gbl] 

a,b=l:bj^a 



Zr = X*"" 



eer /er 



where 



Gt 



Us 

eedf 



o(eJ) 



(7.27) 



(7.28) 
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The function Wf^^{g) is defined in the following fashion: 



WfAa) and Wf^owf^{g)= / dhwf^{gh ^) Wf^{h). 



G 



(7.29) 



In simple terms, not all the face amplitudes turn out to be the same. But there is an obvious 
case when they do all coincide, that is, when Wf^ is stable under group convolution. A nice 
way of ensuring this is to attempt to find the 2d fixed point face amplitude (no matter what 
the dimension of the GFT in question). This is still a highly non-trivial problem and we could 
ask for a slightly more lenient property, say, that the functional form of the wedge amplitude 
is stable under convolution. To take Yang-Mills as an example, the Wilson amplitude does not 
suffice since: 

w;y„((7) =exp(-a9^ex''(5)) (7.30) 

does not maintain the same functional form under convolution of the wedges. But it is well 
known that the heat kernel amplitude is stable under convolution (at the level of each face): 

Wf^ig) = Y,dp x'{g) exp (-a C) =^ Wf,#{g) = Y.d, x'ig) exp (-(#/, x a) C) (7.31) 



peG 



where dp, C^, are the dimension, quadratic Casimir and character of the representation of G 
labelled by p\ a is the coupling constant and is the number of wedges in /. Although the 
weight maintains the same functional form under convolution, one acquires a dependence 
for each face and i^f^ differs for each face. This is fine in the pure gauge theory setting since the 
heat kernel amplitude generically comes with a factor of the area of the face. If we choose all 
wedges to have unit area, then is exactly the area of the face that we wish to pick up. Of 
course, if one wishes to couple this theory to gravity at the GFT level, then one would like this 
area factor to be generated dynamically from the gravity side of the theory. As yet, there are 
no proposals for such a mechanism existing in the literature. Perhaps what is most interesting 
at this stage, is that knowing about the 2d theory can guide us, at least initially, in setting out 
the n-dimensional theory. The action for the gauge model above is: 



dga dga 



a=l 



(f>{9l,---,9n) 



n ^(9ag~a 



a=l 



(t>{9l,---,9n) 



-A 



QXn(n + l) 



n+1 



a,b=l:b^a 



'n+1 



1) 



n+1 

n ^(9a,b9b}) 
a,b=l:b^a 

(7.32) 

where K o Wf^{g) = Wf^ o K[g) = 6{g). There is some freedom here, we could replace the 
5- functions in the vertex operator with Wf^ instead and leave the propagator unchanged, but we 
prefer to act in this way due to its similarity with another class of plaquette theories. These are 
modifications of BF-iheoxy which can be viewed as tentative proposals for higher dimensional 
quantum gravity, the EPRL/FK group field theories. These have a similar functional form to 
()7.32p . obtained by modifying the propagator. Details of this construction can be found in |105| 
and references therein. 



8 Non-local spin foams 

One of the main open problems in quantum gravity is to study the large scale behaviour of the 
models. For this renormalization group methods could be clearly useful. However the application 
of renormalization group methods to spin foam models currently discussed as quantum gravity 
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candidates [311 [32J is hindered not only by many conceptual issues but also by the tremendous 
complicated amplitudes emerging from these models. Here 'toy spin foam models' could help to 
develop new techniques and to investigate the statistical field aspects of these models, some of 
which might be independent of the specifics of the amplitudes. For instance the question whether 
and how to sum over topologies of spin foams might not depend on the chosen group underlying 
the model. Also note that there are a number of quantum gravity approaches which rather work 
with quite simple (vertex) amplitudes [Ml 13 [9] , in particular [9] in which the question of the 
large scale limit can be addressed. The models proposed in this work can be seen as small spin 
cut-off models of the full theory. The hope is that for these models one can probe the many 
particle (and small spin) regime which is in contrast to the very few particle and large spin 
regime considered up to now. 

There has been some very interesting work exploring coarse graining in the spin foam context 
[1211 1122) . Also the related concept of tensor networks |123j . a generalization of spin nets 
introduced here, has been developed as a tool for coarse graining (to specify topological order in 
condensed matter systems) . Mostly these however concentrate on frameworks and truncations in 
which the local form of the spin foams does not change. One might however need to accommodate 
also for non-local^ couplings, in particular if one attempts to regain diffeomorphism symmetry, 
which is broken by the discretizations employed by many quantum gravity models [16l l70] . 

Indeed applying block spin transformations to a non~topological lattice gauge system one 
can expect to obtain a partition function of the form 

oo 

^ = E n n wi,^...iMihh,---,hij (8.1) 

9e M = 11i,...,Im 

where hi. are holonomies around loops (that is around plaquettes but also around other surfaces 
made of several plaquettes) and wi^^,,,ij^j is a class function of its arguments, describing possible 
couplings between (Wilson) loops. A character expansion would not only introduce representa- 
tion labels on the basic plaquettes but also on all the other surfaces encircled by loops appearing 
in (jS.ip . The resulting structure is akin to a two-dimensional generalization of a graph. Graphs 
would appear as effective descriptions for the coarse graining of the Ising like models discussed 
in section [2l Such non-local 'spin graphs' would be generalizations of the spin nets introduced 
there. Similarly, graphs have been introduced in [12] to accommodate non-local couplings and 
to describe phase transitions between geometric (i.e a phase where for instance a space time 
dimension can be defined) and non-geometric phases. We leave the exploration of the ensuing 
structures for future work. 

9 Outlook 

In this work we discussed several concepts and tools which arose in the spin foam, loop quantum 
gravity and group field theory approach to quantum gravity and applied these to finite groups. 
We encountered different classes of theories, one is the well known example of Yang Mills like 
theories another are topological BF theories. The latter are also well known in condensed matter 
and quantum computing. A third class of theories are the constrained models discussed in 
section [5^ which mimic the construction of the gravitational models. These are only applicable 
for non-Abelian groups as for Abelian groups the invariant Hilbert space associated to the edges 
is always one-dimensional and cannot be further restricted. We plan to study these theories in 
more detail in the future, in particular the symmetries and the associated transfer or constraint 
operators. The relative simplicity of these models (compared to the full theory) allows for the 

^^These couplings are in general exponentially suppressed with respect to some non-locality parameter, like 
distance or size of the Wilson loops. In this case one would still speak of a local theory. 
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prospect of a complete classification of the choices for the edge projectors and hence of the 
different constrained models. For the study of the translation symmetries in the non-Abelian 
models it might be fruitful to employ a non-commutative Fourier transform |124l 11161 I117j . 
This would define an alternative dual for the non-Abelian models, in which the edge projectors 
carry delta function factors, however defined on a non-commutative space. 

We also mentioned the possibilities to obtain gravity like models by breaking down the trans- 
lation symmetries of BF theory. This could be particularly promising in the canonical formalism 
where one can be guided by the form of the Hamiltonian constraints for gravity. This strategy is 
also available for the Abelian groups. In particular for Z2 it is in principle possible to parametrize 
all possible Hamiltonians or partition functions and to study the associated symmetry content. 
See also the universality result in [48] . Hence there might be a definitive answer whether gravity 
like models exists, i.e. 4D (Z2) lattice models with translation symmetries based on the 4-cells 
(or the dual vertices). 

Finally we hope that these models can be helpful in order to develop techniques for coarse 
graining and renormalization of spin foam models and in group field theories. Here the connec- 
tion to standard theories could be exploited and techniques can be taken over from the known 
examples and with adjustments be applied to quantum gravity models. To this end the finite 
group models could be an important link and provide a class of toy models on which ideas 
can be tested more easily. For instance the connection between (real space) coarse graining of 
spin foams and renormalization in group field theories, which generate spin foams as Feynman 
diagrams, could be explored more explicitly than in the (divergent) SU (2)-based models. 

It will be particularly interesting to access the many-particle regime, for instance using 
Monte Carlo simulations, which for the full models is yet out of reach. In the ideal case it might 
be possible to explore the phase structure of the constrained theories and to study the symmetry 
content of these different phases. 
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